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We study the critical behavior of a driven interface in a 
medium with random pinning forces by analyzing spatial and 
temporal correlations in a lattice model recently proposed by 
Sneppen [Phys. Rev. Lett. 69, 3539 (1992)]. The static and 
dynamic behavior of the model is related to the properties of 
directed percolation. We show that, due to the interplay of 
local and global growth rules, the usual method of dynami- 
cal scaling has to be modified. We separate the local from 
the global part of the dynamics by defining a train of causal 
growth events, or "avalanche", which can be ascribed a well- 
defined dynamical exponent zi oc = 1 -+- £ c ~ 1.63 where £ c is 
the roughness exponent of the interface. 



I. INTRODUCTION 

The behavior of a driven interface subjected to 
quenched random forces plays an important role in the 
ordering kinetics of impure magnets and other domain 
growth phenomena Q . The driving force F can be real- 
ized by a magnetic field, pressure or chemical potential fa- 
voring the growth of one of the coexisting phases. If F is 
weak, the interface is typically pinned in one of the many 
locally stable configurations. In this case growth is pos- 
sible only through thermally activated hopping, which is 
an extremely slow process at low temperatures. However, 
when F exceeds some critical value F c , all metastable 
states disappear and the interface is then free to move 
even at zero temperature. The depinning of the inter- 
face at F c can be considered as a critical phenomenon 
where characteristic quantities show power-law behavior, 
e.g. the velocity of the interface is expected to scale as 
v ~ (F — F c ) e , where 9 is a critical exponent. It is known 
that charge-density waves pinned by impurities exhibit 
very similar behavior The dynamical behavior asso- 
ciated with the disappearance of metastable states (i.e. 
avalanche) as F is increased towards F c has been a sub- 
ject of recent interest in the study of systems far from 
equilibrium ^|). 

A plausible continuum description of the interface dy- 
namics is given by the following equation with a Kardar- 
Parisi-Zhang (KPZ) nonlinear term, 

!>h — X -(Vhf + F-r,(x,h), (1) 



dt 



where h(x,t) is the height of the interface. Unlike the 
original KPZ-equation for Eden-type growth processes, 
here 77 (x, h) is a quenched random force with short range 
correlations. In the case A = 0, all critical exponents 
describing the depinning transition have been calculated 
recently in a functional renormalization group treatment 
close to four interface dimensions For growth in an 
isotropic medium, it is plausible that the A-term is not 
present when the interface moves with vanishing velocity. 
However, the term can be present for anisotropic growth. 
Most solid-on-solid type models are expected to be in the 
latter category. 

In 1+1 dimensions, there is a particularly simple class 
of lattice models which exhibit a critical depinning tran- 
sition ]9|,[To|. The mechanism for pinning is the directed 
percolation of cells with pinning forces r\ greater than 
the driving force F. The threshold force F c needed to 
depin the interface is then simply related to the critical 
percolation density p c of such cells by F c = 1 — p c - The 
roughness exponent of the pinned interface is equal to 
that of the critical percolation cluster, £ = Q c ~ 0.63. 

In a separate development, Sneppen introduced a sim- 
ple model (model B of Ref. |Ll] , hereafter referred to as 
the Sneppen model) to examine the interplay between 
local and global rules of growth in determining interface 
roughening and temporal correlations. He found numer- 
ically in (l+l)-dimensions, that the roughness of the in- 
terface also obeys the scaling of a string on a critical 
directed percolation cluster. 

As we explain below, the same roughness exponent for 
the two types of models is due to the fact that the same 
geometrical constraint (the "Kim-Kosterlitz" condition) 
is invoked in defining a locally stable configuration p2[ . 
There are however differences in the way the interface 
is driven. Although the pinning and elastic forces on a 
given site in the Sneppen model are local as in eq.(l), each 
time the site on the interface with minimal ij — 7? m i n is 
selected globally and made to grow by one unit in height. 
Neighboring interface sites then adjust themselves to re- 
cover the geometrical constraint. In the language of a 
uniformly driven interface, such a rule corresponds to in- 
creasing F just above ry m i n to make one site unstable, and 
then quickly set F to a much smaller value to prevent an 
avalanche taking place. 

The Sneppen rule allows one to sample a particular 
sequence of interface configurations, each of them being 
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a metastable configuration at some value of F. Succes- 
sive configurations in the sequence differ only by an in- 
finitesimal amount, i.e. a few sites (about 4) which have 
moved when one site is made unstable. The situation 
here resembles that of an interface at a finite tempera- 
ture and driven by a uniform force F far below F c : Local 
irreversible motions are made possible by thermal fluc- 
tuations, but the interface stays always close to some 
metastable state. 

Further numerical studies of the model by Sneppen and 
Jensen ]l3| revealed interesting spatial correlation of suc- 
cessive growth events. In addition, they have found that, 
in the saturated regime, the height advance at a given 
column exhibits complicated scaling with time, with ex- 
ponents varying with the moment considered. 

The aim of the present paper is to analyze the spatial 
and temporal correlations in the Sneppen model and try 
to relate the observed scaling behavior to the properties 
of directed percolation. It turns out that some of these 
correlations depend on the value of 7/ m i n (or equivalently 
F) at a given moment. Since rj m i n fluctuates in time, for 
these correlations the temporal translational invariance is 
lost. This is an important feature of the dynamics based 
on global rules. 

Another significant consequence of the global rules is 
that growth is no longer homogeneous in space. At a 
given moment, only a small part of the interface is mov- 
ing. The growth events that follow may either closeby, or 
far away. Due to this property, the usual method to per- 
form dynamical scaling, based on the presumption of a 
single growing correlation length, should be modified. In 
this connection we found it useful to distinguish growth 
events which are closeby and hence bear strong correla- 
tions, from those which are far apart. We observe that 
a train of growths events, started with some rj 1 ^^, prop- 
agates laterally with a well-defined dynamical exponent 
zioc, which can be related to the roughness of the di- 
rected percolation cluster: zi oc = 1 + £ c ~ 1.63. In the 
context of a driven interface, this motion can be thought 
of as an avalanche. We found that the distribution of 
avalanche sizes obeys a power-law decay up to a size re- 
lated to r]® nin . The spatial-temporal correlations between 
successive growth events, on the other hand include both 
local and global motion. 

The paper is organized as follows. In section II we re- 
call the definition of the Sneppen model and present a 
theorem which relates the stable (static) configurations 
of the Sneppen model to directed percolating strings. In 
Section III we define the avalanches (causal events) and 
determine their dynamical behavior as well as their size 
distribution. In Section IV the spatial-temporal correla- 
tions are investigated by considering the distribution of 
lateral distances between successive growth events. Sec- 
tion V contains conclusions and a summary. 



II. DISTRIBUTION OF PINNING FORCES AND 
ROUGHNESS 

We first review the definition of the Sneppen model 
[0|. Each cell (i,h) on a square lattice is assigned a 
random pinning force r)(i, h) uniformly distributed in the 
interval [0, 1). The interface is specified by a set of inte- 
ger column heights hi (i = 1, —,L) with the local slope 
constraint \hi — hi-\\ < 1 for all hi ("Kim-Kostcrlitz" 
condition). Growth hj — ► hj + 1 proceeds at the site j 
where the pinning force rj(j, hj) = r/ m j n is the minimum 
among all interface sites, followed by necessary adjust- 
ments at neighboring sites until the slope constraint is 
recovered. The growth rules are illustrated in Fig.l. 

Since we want to relate the behavior of the interface in 
the Sneppen model to directed percolation, we first recall 
some of the properties of directed percolation clusters 
0. When the density p of occupied sites is less than 
some threshold value p c , a typical cluster of occupied 
sites connected horizontally or diagonally extends over a 
distance of the order of £|| in the parallel direction and a 
distance of the order of £j_ in the perpendicular direction. 
For p > p c , there appears a directed percolating cluster 
which extends over the whole system. This cluster has 
a network structure of nodes and compartments, where 
each compartment has an anisotropic shape similar to 
the connected clusters below p c , characterized by £|| and 
£_L- On both sides of the percolation transition, the two 
lengths have the power-law behavior 

su-iP-Pcr 11 , tL~\p-pc\- uj -. (2) 

Series calculations give v\\ = 1.733±0.001 , v± = 1.097± 
0.001 HI and p c = 0.5387 ± 0.003 @. The roughness 

of a percolating string scales as £j_ ~ £y X j i-e. the 
roughness exponent £ c = v±/vh ~ 0.63. 

In Ref. Jl2j we proposed to study the distribution P v {rj) 
of pinning forces rj(i,hi) at the interface and the probabil- 
ity distribution P m (i]min) during growth. When starting 
at t = with a flat interface hi = 0, the forces hi) 
are equally distributed. During the transient regime the 
interface becomes rough and since always the smallest 
pinning force ?y m in is updated, the sites with small r}{i, hi) 
get rare. This in turn implies that the typical value of 
the selected 7y m i n increases with time. The distributions 
P p {rj) and -P m (f?min) shown in Fig. 2a were recorded in 
the transient regime in the time interval L/A < t < L/2 
for a system of size L = 8192. The peak of P m (^min) 
moves to the right with increasing time and thereby " eat- 
ing up the store" of small rj(i, hi) which were present at 
t = for the flat interface. As long as there were no 77 m ; n 
larger than a value rj u (t), the distribution P p (ji) is still 
a constant for r\ > r/ u (t). (In Fig. 2a r) u (t) w 0.3.) In 
the next paragraph we show that in the thermodynamic 
limit, there will never be an 7? m i n larger than a critical 
value F c — 1 — p c ~ 0.461 and therefore the transient 
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regime ends when ?y m i n first comes close to F c . When 
the peak of P m (Vmin) approaches F c , its height vanishes 
and the distribution becomes stationary, which we show 
in Fig. 2b together with P p {rj) in the saturated regime. 
Since ?7 m in < F c , the stationary distribution P p [rj) is flat 
for rj > rj u = F c in the limit L — ► oo. 

To see that the growing interface always has r7 m m < 
F c in the thermodynamic limit we first note that every 
interface configuration satisfying the slope constraint is 
a path on a directed percolation cluster of sites with rj 
greater than or equal to ?y m i n . Such a path only exists if 
the density 1 — ry m i n of these cells on the lattice is greater 
than the critical percolation density p c , i.e. r/mm < 1 — pc 
for all interfaces. Paths on the infinite critical percolating 
cluster have the largest T] m i n = F c = 1 — p c ~ 0.461. 

In a numerical simulation with a finite system however, 
we see only a part of an infinite critical percolating clus- 
ter. Thus, an interface which traces out a critical path 
can have a value ?y m i n slightly larger than F c and the dis- 
tribution P m (?7 m in) in the saturated regime in Fig. 2b is 
not exactly zero for ?7 m i n > F c . The motion of the flat 
interface at t = to the first critical path corresponds to 
the transient regime (see Fig. 2a). 

In the following we show that an interface in the tran- 
sient regime as well as an interface which already crossed 
a critical path, is driven to configurations with succes- 
sively increasing 77 m i n and thereby approaching the next 
path on a critical percolating cluster. We first introduce 
a few notations. Assuming that the random forces rj 
are real numbers, there will not be two interface con- 
figurations sharing the same 7y m i n because an interface 
is always updated at the site with rj m i n . Thus each in- 
terface configuration {hi} can be characterized by its 
7y m i n and is denoted by H[rj m i n ]. An order relation is 
defined by H[r)£ in ] > H[rjg in } if hf > hf for all i 
and #[?7miiJ H[ 7 lniin]- Consider an interface fffamin] 
at a time t = to and choose a real number c with 
< in <c<F c . 

We next show that the growing interface at times t > to 
has to overlap completely with the closest percolating 
path which has rfaiin > c. This closest path is defined 
by H[r]f nin ] = min{ fffomin], such that r; min > c and 
fffamin] > ^[^min] }• There can be many percolating 
pathes which share the same but only the lowest 
path H [77^] will be realized by the interface as can be 
seen as follows. 

The growing interface configurations H [r7 m ; n ] at times 
t > t which are below H[r]^ in ], H[rf^ in ] < H[rj min ] < 
^famiJ) a11 nave an Vmin < Vmm b y tne definition of 
^[Vminl- (This is also true for interfaces H[rj nlin ] which 
already overlap in part with H [t^J , no matter if the site 
with 77^ in is already on the interface H[r} m i n ].) Therefore 
these r/ m i n will all be updated before 77^ in and hence the 
interface will grow until H[r/ mul ] — iJfjy^J and then up- 
date the site with r7^ in because any interface site can 



advance at most by one unit in height upon each updat- 
ing. 

The motion with c = rj® lln we will call later an 
"avalanche" (see Sec. III). If we choose however c = F c , 
it follows that all pathes which are the minimum on a 
critical percolation cluster, i.e. -ff^^m] — H[F C ) in the 
thermodynamic limit, act as "checking points" where the 
interface has to go through. There, a "snapshot" of the 
distribution P p (rj) is exactly a step function. It can be 
seen from Fig. 2b that the intermediate interface config- 
urations between two checking points are also close to 
critical, in the sense that only a small percentage (which 
we expect to vanish in the thermodynamic limit) of sites 
on the interface have rj < F c . Hence the interface in 
the Sneppen model has approximately the roughness ex- 
ponent C of the critical directed percolating path with 
C ~ C - 0.63. 

We have seen that the interface H [r] m i n ] is successively 
driven to configurations iTfi^iJ with rj^ lhl > 77 m j n and 
thereby eliminating r}(i,hi) < f^in 011 ^ ne interface. This 
process (see Fig. 2) may be perceived as "self-organizing" 
part of the approach to criticality, where T)mm = f° r in- 
finite systems. In comparison, the interface in the model 
of Ref. is pinned by the critical percolating cluster 
when the driving force F is tuned to its threshold value 
F c = 1 - Pc . 

The roughness exponent £ can be measured by the 
equal time height- height correlation function C(r) — 
([h(r + r' ,t) — h(r' ,i)] 2 } ~ r 2 ^ where the overbar and 
the angular brackets denote the spatial and the con- 
figurational average, respectively. For a system of size 
L = 8192 we found a roughness exponent ( = 0.655 ± 
0.005 which is somewhat larger than the critical value 
0.63. However, the measured exponent £ varies sys- 
tematically with the system size: For L — 900 we get 
C = 0.665 ± 0.005, whereas for L — 65536 we found 
C = 0.648 ± 0.005. An explanation is that in a finite 
system the distribution P p {rj) is not exactly a step func- 
tion and the motion between two critical clusters yields 
an effective exponent larger than 0.63 similar to a moving 
interface of Ref. |1. In our simulations, P p (ri) approaches 
the step function with increasing system size. Thus we 
expect that the measured discrepancy to the exponent of 
the percolating cluster £ c ~ 0.63 vanishes in the thermo- 
dynamic limit. This is also supported by a simulation 
where we measured C(r) only when ?7 m i n is close to F c . 
In this case £ = 0.64 ± 0.01 which is consistent with the 
expected critical value. 

III. CAUSAL EVENTS: AVALANCHES 

After a transient regime, when rj m { n first comes close to 
F c , the interface in the Sneppen model exhibits a steady- 
state critical behavior (saturated regime) , which allows a 
convenient study of the dynamics at criticality. However, 
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the behavior is complicated, caused by the interplay be- 
tween the local adjustments due to the slope constraint 
and the rule that the growth site with rj = ?7 m i n is cho- 
sen among all interface sites. Successive growth events 
can be far apart and the motion is therefore inhomogc- 
neous in space. Hence, a single growing correlation length 
does not exist and the usual dynamical scaling has to be 
considered with care. The realized sequence of growth 
sites after a time t = to depends on the globally chosen 
value r? m in(t = to), which is responsible for the growth 
inhomogeneity in space. Thus we will separate the local 
part from the global part of the dynamics by defining an 
" avalanche" , which has the property that the sequence 
of growth sites inside the avalanche does not depend on 

??min(i = t ). 

An avalanche is defined by a sequence of growth events 
(including the necessary adjustments due to the slope 
constraint) started at any integer time t = to where 
??min(i =to) is denoted by ?7^ in . This avalanche is termi- 
nated at the first time r when the successive n m i n (T + 1) 
is larger than ?7^ in , i.e. for all times t with t < t < t 
the growth events have w m i n (t) < fy^m- We call this 
causal events because the avalanche consists of a train 
of growths events which are all induced by local adjust- 
ments in hi due to the slope constraint after t — to- 

To see in what sense the sequence of growth sites inside 
an avalanche is independent of ??^ in , we consider at t = 
to two identical interface configurations A and B with 
the same random forces r)(i, h) above the interface but 
with different r](i, hi) at the interface such that r/m in [^4] < 
^minl-E?] at the same site j. For the configuration A there 
exist forces r](i, hi) with r^in^ 4 ] <vih h i) < Vmini B ] but 
for interface B all r)(i,hi) > Vmini^]- Since the random 
forces above the interface are assumed to be the same, 
all Vmin of the growing interfaces A and B inside both 
avalanches are identical because they are all induced by 
identical local adjustments after t = to- When however 
at time t a + 1, < in [A] < n min {T A + 1) < J?° in [-B], i.e. 
when avalanche A is terminated, this new growth site of 
interface A can be far away, but for the interface B it has 
to be still induced by the local adjustments of avalanche 
B because there were no r](i,hi) < ^^[B] at t = to- 
We have seen that although the random environment at 
(and below) the two interfaces A and B is different, the 
motion inside the avalanche is only influenced by the local 
adjustments after t = to- 

A size s of the avalanche can be defined by the num- 
ber of growth events including the necessary adjustments, 
s = J2i(hi( T ) — hi(t )), and awidth by max{ i, such that 
hi(r) > hi(t ) } — min{ i, such that hi(r) > hi(to) }■ 

From Sec. II we know that an interface with T)^ in is 
driven to configurations with r^ in > ?7^ in . At t = to 
with n m i n (t = t ) = ?7min> a P ar t of the interface starts to 
move which was "pinned" just before to by a path with 
i](i,hi) > ^ in . This part of the interface moves through 
a compartment of the percolation cluster to the next path 



which again pins the interface with ij(i,hi) > c = rj^ in 
(see Sec. II). Since the compartment has a height of the 
order of £±{Vmin) an( l a width of the order of £||(?7min)' 
it is a natural conjecture that the width of the avalanche 
scales with ^|i(^in) ~ Wc — ^mml ^" an d the size is at 

mOSt £||« ;n ) * a« in ) ~ [F c - vLn]- {U±+ ^- 

Note that at every (integer) time an avalanche is 
started and big avalanches can contain smaller ones. 
Thus successive growth events inside a big avalanche can 
be quite far apart (jumps between small avalanches), but 
all events are inside the correlation length ClK^min) ~ 
[F c — ^mm]^ 11 • I n this sense an avalanche is "localized" 
and we call the corresponding motion " local dynamics" . 
The presumption of dynamical scaling, that there is only 
a single growing correlation length, is reestablished for 
the local dynamics and we can ascribe a well-defined lo- 
cal dynamical exponent z; oc to the lateral propagation of 
growth inside the avalanche: r <~ ^ ,oc . Since the size of 
an avalanche is proportional to the time r, one has from 

z loc ~ 1 + i/ ± /i/|| = 1 + Cc 1-63. (3) 

In a simulation, this local dynamical exponent can be 
detected by considering the infinite moment of the height- 
height time correlation function in the saturated regime 



C q (t) = ([hi{t + f) - hi{f) - hi(t + 1') - hi(t')}9 J , 

(4) 

which becomes for q — > oo 

Coo(t) ~ {maxi{hi(t + t') - hi(t')}) . (5) 

Since only the column with maximum height-advance 
A/i max contributes to the infinite moment Coo(i), most 
growth events during the time t are involved in the mo- 
tion through a compartment of a percolation cluster of 
height A/i max , i.e. Ah max - £ ± - ~ t 1 " 1 ^- = 
t CAi.c. Thus Coo(t) - &<*> with 13^ = C/z loc = (3 loc . 
Sneppen and Jensen observed a scaling of Coo(t) with 
Poo = 0.41 ± 0.02 which is close to the exponent (3i oc = 
C/zioc — v±/(v± +v\\) — 0.39. Our own simulations give 
/3oo = 0.40 ±0. Of which is in perfect agreement with (3i oc 
if we insert the measured £ ~ 0.655. 

The distribution Ph(Ah, t) of height advances Ah(t) = 
hi(t' + t) - hi(t') is shown in Fig.3a for 2 6 < t < 2 15 
and Ah > 0. Ph has a large peak at Ah = which is 
not shown, i.e. most of the columns i have not grown 
(over 99 percent for t = 2 6 and 70 percent for t = 2 15 
with L = 8192). Next we show that the distribution 
Ph(Ah, t) for the moving columns can be roughly brought 
to a "local" scaling form when Ah is scaled by i /3loc . 
To normalize Ph(Ah > 0,t) we note that the portion 
n m ov of columns which have moved, scales as the correla- 
tion length divided by the system size, n mov ~ t 1 / z,oc /_L. 



4 



Thus one has 



P h (Ah>0,t) 



t l/z 



t 



1-2/3,, 



Ah 



Ah 



(6) 



where T(y) is a scaling function (see Fig. 3b). For fast 
growing columns with large Ah/t^ 1 "" the above argument 
for Ah max applies and the data collapse in Fig. 3b is per- 
fect. For smaller Ah/t 13 '"", however, there are significant 
deviations from the "local" scaling form. 

For the second moment C2 (t) Sneppen and Jensen 
observed a scaling with an exponent 02 — 0.69 ±0.02 [ p^[ . 
However, due to the inhomogeneity in growth, the appli- 
cation of dynamical scaling is questionable as mentioned 
above. The deviation of the effective exponent 02 from 
the scaling of the local dynamics is caused by the fact 
that only a small part of the interface (n mov ~ t 1 ^ 10 " /L) 
has moved for short times t. The observed value for 02 
can be explained by using that Ah scales roughly with 
fPioc f or Aft, > (see Fig. 3b). For general integer q we 
write 

7I/9 



C q (t) = ( (Ah -Ah) 



Ahi - qAhi- 1 Ah 



Ah 



1/'/ 




q l t (9-l)ft..l 



(Ah(t) = t/L). In the observed scaling regime t <C L Zloc 
and therefore the main contribution comes from the first 
term. Thus we have 



C q (t) ~ t l3loc+1 / qz '° c n*> t 1 +( 1 -q)/ c i z i°c 



(7) 



i.e 



q ~ 1 + (1 - q)/q(l + C) and 2 - 0.70 which agrees 
with the simulations, while for q — > 00 , q — > 0\ oc . 

Next the avalanche size is investigated. We observe 
that the distribution P av (s, rp^^) of the avalanche size s 
for a given T]^ ain shows a power-law decay with an expo- 
nent k~ 1.25 ±0.05 up to a size s ~ C|| (^Sun) *^(Vmin) 
and then drops to zero for s > so (see Fig. 4a). Thus the 
avalanche size distribution obeys the scaling form 



Pav(s) 



so 



(8) 



with $(y) = const, for y < 1 and a rapid decay for 
y > 1. The good data collapse onto the scaling form 
eq.(8) in Fig. 4b supports our picture that an avalanche 
corresponds to the motion through a compartment of a 
percolation cluster, which is characterized by the lengths 
£1 and £|| . 



IV. SPATIAL-TEMPORAL CORRELATIONS 

In this section we try to understand the spatial- 
temporal correlations between successive growth events. 
To this end we investigate the probability distribution 
P co (x,At), where x is the distance parallel to the inter- 
face between growth events which occur after a time At. 
Sneppen and Jensen |l3| observed that P co {x, At) is con- 
stant for sufficiently small x and a has power-law decay 
with an exponent 7 = 2.25 ± 0.05 above a value x c which 
increases with At (Fig. 5a). 

We next explain that the behavior P co (x, At) = const 
corresponds to the dynamics of causal growth events in- 
side an avalanche. The local adjustments due to the 
slope constraint after the avalanche has started induce 
randomly distributed rj(i,hi). During the avalanche all 
f?min(i) < Vmin are taken from these newly appeared 
r)(i,hi). Thus, the rj m i n (t) are randomly distributed 
in space, i.e. the distance between successive growth 
events is also equally distributed as long as At is smaller 
than the duration of the avalanche r, i.e. as long as 
x < CllC^min)- Therefore we can cast P co (x,At) into the 
scaling form 



Pm 



(9) 



with zi oc = 1 + £ and ~^(y) — const, for y < 1 and 
$ y~T f or y > 1 The scaling form eq.(9) with a 
satisfactorily data collapse is shown in Fig. 5b. 

We see that the spatial-temporal correlations depend 
on the value rj m i n . For 77 m j n close to F c , x is equally 
distributed even for large x. For small rfaiin on the other 
hand, P co (x 1 At) obeys a power-law decay also for small 
x, i.e. it is more probable that successive growth events 
are closeby. Thus, for these correlations the temporal 
translational invariance is destroyed. 

The value of the exponent 7 can be understood by 
considering the conditional distribution P co for fixed ?7 m i n 
of the first event, which we denote by P co (r) m i a ,x, At). 
We express P co (x, At) by 

P co (x, At) = J dr) min Pa, (r) min ,x, At) P n (ry roin ). (10) 

From our simulation (see Fig. 2b) we see that the 
probability distribution P^(?/ m i n ) ~ (F c — ?/ m i n ) close 
to F c . From the above discussion we know that 
Pco(Vmin,x,At) ~ 1/£|| if £|| > x. Thus the integral 
eq.(10) takes the form 

P co (x, At) = [ di] min -!- (F c -r] min ) 

from which one obtains P co {x, At) ~ x" 1 with 

7 = 1 + (2/i/||) ^ 2.16. (11) 
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This is quite close to our numerical result 7 = 2.20 ± 
0.05. We have also directly measured the distribution 
P co (f?min, x, At) to check our assumptions. We found that 
Pco (f?min, x, At = 1) first decreases for small x due to 
a high probability for choosing the next ?7 m i n from the 
newly appeared r\ from the local adjustments. However, 
to explain the exponent 7 we are interested in large x 
(and thus in large £11), for which we indeed observe a 
constant P C o(r)min, x, At) for x < £|| (f7 min ). 
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V. CONCLUSIONS AND SUMMARY 



We have analyzed a number of spatial and temporal 
correlations in the Sneppen model jll 13 1. The imposed 
slope constraint \hi — < 1 is the reason for the 

relation of the static and dynamic behavior to the prop- 
erties of directed percolation. The roughness exponent 
of the interface in the Sneppen model and of the pinned 
interface in the model of Ref. || is equal to that of a 
percolating string, C c ~ 0.63. 

The difference between the two models is, however, 
that in the model of Ref. || the interface is driven by a 
uniform force whereas in the Sneppen model there is a 
sclf-tuned driving force which keeps the interface at the 
onset of steady-state motion and therefore the interface 
shows critical behavior. This is achieved by the rule, that 
the site grows, which has the weekest pinning force r) m i n 
among all sites of the interface. This induces a nonlocal 
part in dynamics. As a consequence, the motion of the in- 
terface is inhomogeneous in space and the methods of dy- 
namical scaling are not applicable in a direct way because 
there is no single growing correlation length. Thus we 
have separated the local from the global part of the mo- 
tion by introducing an " avalanche" , and assigned a well- 
defined dynamical exponent z\ oc = 1 + £ c to the lateral 
propagation of the growth inside an avalanche. We found 
that the size distribution of the avalanches started with 
r /min nas a power-law decay with an exponent k ~ 1.25 
up to a size £||« in ) * UivldJ- 

The spatial-temporal correlations where investigated 
by the probability distribution P co (x, At) which shows a 
crossover from a behavior determined by causal growth 
events (P co (x, At) — const.) to a power-law decay with 
an exponent 7 ~ 2.2 which can be also related to expo- 
nents of directed percolation. We have seen that for the 
distribution P co the temporal translational invariance is 
lost, which is due to the global part of the dynamics. 

Upon completion of the paper we became aware of 
an independent work by Z. Olami, I. Procaccia, and R. 
Zeitak where ideas similar to ours have been developed. 
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FIG. 5. (a) Logarithmically binned distribution of dis- 
tances x between two growth events, where P co (x,At) is the 
density of x in the range [x, 2x). The two growth events occur 
after a time At (1 < At < 256, L = 8192). Curves with a 
larger At have a wider plateau. The power-law decay has an 
exponent 7 = 2.20 ± 0.05. (b) Scaling plot of the distribution 
P co (x,At) eq.(9), where x is scaled by At 1 ^ 1 ^' using the 
measured value £ = 0.655. 



FIGURE CAPTIONS 



FIG. 1. An interface configuration (filled circles and full 
lines) is updated (big circles) at the site which has the smallest 
pinning force ?7 m in. Four possible growth events are shown 
including local adjustments (small circles and dashed lines) 
due to the slope constraint \hi — hi-i\ < 1. Time is measured 
in units of such growth events. Each interface configuration 
may be considered as a directed percolating path. 



FIG. 2. Distribution P P (rf) of all pinning forces r](i, hi) at 
the interface (full line) and probability distribution P m (?7 m m) 
(dashed line) for a system of size L — 8192. (a) The distri- 
butions in the transient regime are averaged over the time 
interval L/4 < t < L/2 and over 1000 independent runs, 
(b) P P (rj) and P m (»?min) averaged over time in the saturated 
regime. 



FIG. 3. (a) Distribution Ph(Ah,t) of height advances 
Ah > for time differences 2 6 < t < 2 15 (L = 8192). Higher 
curves correspond to larger t. (The same plotting symbol is 
used for data at a given time t.) In addition, there is a large 
peak of Ph at Ah — which is not shown. The curves are 
normalized if the peak at Ah = is included, (b) Scaling plot 
eq.(6) with the local exponent /3; oc ~ 0.40. 



FIG. 4. (a) Logarithmically binned distribution of 
avalanche sizes s, where P av (s, r; m i n ) is the density of events 
in the range \s,2s). Lower curves correspond to smaller 
r? mi „. (0.30 < J7 min < 0.46, L = 8192). P av (s, Vmia ) has a 
power-law decay with an exponent k ~ 1.25 ± 0.05 up to a 
size so ~ £|| (??min) *Ci( 7 7min) an( l then drops to zero for s > so- 
(b) Scaling plot of the avalanche size distribution (eq.8). 
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